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Оптимальный алгоритм для решения 
интегральных уравнений Гаммерштейна 
в системах компьютерной алгебры 


В работе построен алгебраический алгоритм для вычисления алгебраического многочлена у, порядка 
п = М в системах компьютерной алгебры. Этот многочлен аппроксимирует решение у(х), х Е [а,6] 
интегрального уравнения Гаммерштейна. Эта аппроксимация оптимальна в пространстве Са. 


Актуальность темы исследования. Обыкновенные дифференциальные и инте- 
гральные уравнения являются классическим аппаратом математического моделирова- 
ния [1]. Уравнения Гаммерштейна вида 

У] + = Н»], Г, Н : раупопиа! —> рой1упопиа|, # е роупопиа1 (1) 
используют в математических моделях более часто, чем интегральные уравнения дру- 
гих типов. Маре, Маетайса, Ма сад, МаНаб, АР$ и другие компьютерные системы 
стали естественной средой математического моделирования. Эти системы выполняют 
символьные (аналитические) преобразования композиции специальных математических 
функций и решают обыкновенные дифференциальные уравнения, являются система- 
ми компьютерной алгебры (СКА). 

Задача работы. Построить алгебраический алгоритм для вычисления в СКА 
алгебраического многочлена у, порядка и Е №. Этот многочлен аппроксимирует 
решение у = у(х), х Е [а,в] интегрального уравнения Гаммерштейна вида (1). Эта 
аппроксимация оптимальна в пространстве Сиьы, в этом пространстве ограничен 
коэффициент оптимальности алгоритма 


Си( авогийт, (1), Сань ) = [у — У [аьы / Ш о,...е_п | У) — (со +...+ си Х') аа. 
Актуальность задачи. СКА не решают интегральные уравнения вида (1). 


1. Алгоритм 1 


Вход. Уравнение (1) — 1а5^ (1), отрезок [а,6], параметр п. 

Преобразования. 

1. Вычислить начальное приближение к решению уравнения (1) у=0. 

2. Для 5 = 1,2,... вычислить: 
— преобразование / = /(х, ун;-1 ) многочлена у»„,;.! функцией Х (1), 
— алгебраический многочлен, интерполирующий функцию Х, в узлах Чебышева на от- 
резке [а,6] 

Е; = Ч ]=Ч п[{ 5(а+(Ь-а) (1- со п/п ))/2), 1=0,...м }], (2) 

— преобразование многочлена (2) интегральным оператором Н[у] уравнения (1), 
— алгебраический многочлен 


& =Н БВ], (3) 
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— линейное интегральное уравнение с многочленными коэффициентами (ЛИУМК) 
МУ +8 =8» (4) 
— решение ЛИУМК (4) на отрезке [а,6] по алгоритму [2] с параметром п — алгебраи- 
ческий многочлен 
у»; = авогийт_[2] ( Цу| +8 =8., [а,6], п). (5) 
Вычислить предел последовательности алгебраических многочленов (5). 
Выход. Искомая аппроксимация решения у уравнения (1) 


у» = авогийт_1 (ва (1) ‚ [а,Б] , п) = со + си х +... си х" = ШП, о Уи». (6) 
Замечание 1. Алгоритм 1 с итерационным методом Ньютона или методом более 
высокого порядка преобразуют уравнение (1) в последовательность ЛИУМК вида 
НУ+8=Н Ефв,ун-1) У] + СС Н/У &уть- ), $ =1,2,.... 
Решать их по алгоритму [2], как правило, более сложно, чем ЛИУМК (4). 


2. Алгоритм 1 и метод квазилинеаризации на задаче 
У'=ехро›, У(0)=0, ур =0. 


Эта краевая задача является одномерным аналогом уравнения гидродинамики 
Ух + у"и = ехр(у). 
Она является модельной для метода [1] и эквивалентна интегральному уравнению 
у= КГехр() |, К = | (1-х гуф & + а-дх у 4. (7) 
Уравнение (7) и эта краевая задача имеют единственное решение — аналитическую в 
окрестности отрезка [0,1] комплексной плоскости функцию 
у=- \1(2) + 2 Ма(с зес(с (х-1/2)/2), с = зоУе( за (2) = с зес( с/4 )). (8) 
Отличные от нуля коэффициенты Фурье — Чебышева функции у (8) на отрезке [0,1] 
только четные. С возрастанием параметра п = 2? они монотонно убывают к нулю а2 Ку, 
[0,1]) = о(а'), 4 < 0.01 и принимают следующие значения: 
{ а Ку, [0,1] 4-0” = { -0.057, 0.057, 0.00027, 2.1 10°, 1.8 10°, 
1.66 10°, 1.6 107, 1.6 107“, 1.6 1075, 2.1 1078 3. (9) 
Следовательно, для величины наилучшего приближения функции у (8) алгебраичес- 
кими многочленами порядка и в пространстве Сю! справедливо тождество 
Ш 0... || УС®) — (со +... сих’) [аю.л = (1 + 0(1)) [ а2 из] +2(7 ‚ [0,1] |. — (0) 
Вычислительный эксперимент с алгоритмом 1. С возрастанием параметра п 
алгоритма 1, норма погрешности решения уравнения (7) по алгоритму 1 в простран- 
стве Сю] монотонно убывает к нулю, удовлетворяет тождества 
{Пу ур; [< =г = { 0.00034, 3.1 10%, 2 10%, 1.8 1075, 9.7 107°} — (11) 
и скорость убывания тождественна скорости убывания четных коэффициентов Фурье — 
Чебышева (9) функции у (8) на отрезке [0,1] имеют место тождества 
ПУ [0.1 / | 42 1+2, [0,1]) | =1+ ор ирь {ры = { 0.3, 0.5, 0.1, 0.07 }. 
Коэффициент оптимальности метода решения задачи С, ( шео4, 1азК, Са, ) > 1. 
Поэтому из этих тождеств и тождества (10) можно сделать следующее заключение. 
Вывод 1. Коэффициент оптимальности решения по алгоритму 1 уравнения 
(7) в пространстве Сью, асимптотически минимальный 
С‚(авогийт_1, 1а5К (7), Сю, п) = (1 + о(1)) (1+ 092 тр). (12) 
Коэффициент оптимальности метода квазилинеаризации [1]. По методу [1] 
решают (не линейные) функциональные уравнения (145%) на отрезке [а,6] и вычисляют 
сеточную функцию с равноотстоящими узлами 
У = Ув-ф-ауч(хь Фа5К ) }:-0' = тео [1] (ваз , [а,6], 4) = {ус} -0", 
Х0о=@ м=атй,. м6, й= (6 -а)щ. 
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Композиция метода квазилинеаризации и интерполяции алгебраическим многочленом 
сеточной функции у» (вычисленной по методу квазилинеаризации) является методом 
вычисления алгебраического многочлена порядка 4 

Ра у» ] = Ра тефоа [1](а5к, [а,Б], 4)]. (13) 
Этот многочлен аппроксимирует решение у функционального уравнения {45 на отрез- 
ке [а,6]. Норма погрешности метода (13) в пространстве Си] не меньше соответст- 
вующей сеточной нормы. Сеточная норма погрешности метода (13) тождественна 
сеточной норме погрешности метода квазилинеаризации 


пах, с [5] [У — Рау» | > Шахн....4| (7 — Рен) к =х+ | = пах:-о,...4| УС — ущхь 1а5к)|. 
Следовательно, коэффициент оптимальности метода (13) в пространстве С = Сиы не 
меньше 


СеРч[ун, 195, С) > ттах:-0...4| (0%) — ню | / ео... е а | УСО) — (со... сах) [[с. 
Сравнение алгоритма 1 с методом квазилинеаризации [1]. На краевой зада- 
че, эквивалентной уравнению (7), погрешность метода квазилинеаризации [1, с. 42] 


а 


Из этого тождества и тождеств (9) - (12) можно сделать следующее заключение. 

Вывод 2. Коэффициент оптимальности метода (13) на краевой задаче, экви- 
валентной уравнению (7), в 100 000 раз больше коэффициента оптимальности алго- 
ритма [1 на уравнении (7) 

Со(Ри[уь =мо|, ва5® (7), Сю.) > 3 107/ 1.6 107? > 100 000. 

Сравнение алгоритма 1 и методов с насыщением. Метод квазилинеаризации 
и другие сеточные методы решения функциональных уравнений являются методами 
с насыщением. Скорость убывания, с ростом числа узлов, погрешности решения по 
методу с насыщением функционального уравнения с решением аналитической функ- 
цией полиномиальная 

тах;-0,...п| У(ХА — Ув--ауи(х) |= О(и?), р = р(тефо9) = СопзЕ. 
Скорость убывания (с возрастанием параметра п) величины наилучшего прибли- 
жения функции у (8) алгебраическими многочленами порядка п в пространстве Суо,11 
оценивает тождество (10). Поэтому коэффициент оптимальности метода (13) и дру- 
гих методов интерполяции сеточного решения у, многочленами возрастает с ростом 
числа узлов сетки 4 и скорость роста 
СР], 1А5К , Сю, п) = О0(2'), 20, 

Коэффициент оптимальности метода ряда Тейлора имеет такую же скорость роста. 
Из этих тождеств и вывода 1 можно сделать следующее заключение. 

Вывод 3. Согласно критерию, минимум коэффициента оптимальности, на 
краевой задаче, эквивалентной уравнению (7), композиция алгоритма преобразования 
этой задачи в интегральное уравнение и алгоритма 1 предпочтительнее метода (13) 
и метода ряда Тейлора. 


3. Алгоритм 2 — метод Галеркина с возмущением 


Вход. Уравнение (1), отрезок [а,6], параметр и. 


Преобразования. 
1. Вычислить аппроксимацию уравнения (1) по методу Галеркина 
5 Ы у» е Ниа,6]] ] = 5 НЕО Ду» = Ныа,6]) ] ] + 8 ], (14) 


где оператор проектирования 5„ вычисляет частную сумму порядка и ряда Фурье — 
Чебышева функции у на отрезке [а,6] — 5»: Г2(а,б;р) —> Ни[а,6] — с возмущением 
интерполированием (/,[/] (2) по узлам Чебышева. 
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2. Вычислить начальное приближение к решению уравнения (14) у =0. 
3. Решить уравнение (14) по методу простой итерации для 5 = 1,2,...: 
— вычислить преобразование многочлена у; функцией / (1) - } = Хх, уи5-1), 
— вычислить многочлены вида (2), (3) — Е; = Ц} |, ®;, = Ы Е, |, 
— вычислить ЛИУМК вида (4) 


Цу]+8=85, (15) 
— вычислить решение ЛИУМК (15) на отрезке [а,6] по методу Галеркина (14) — ал- 
гебраический многочлен 
Уи; = $0У6( 5,[ [у ЕР, |+82-—8, | =0), 51 : [2(а,6; р) —> На, 6]. (16) 
4. Вычислить предел последовательности многочленов (16). 
Выход. Искомая аппроксимация решения уравнения (1) у, = Шт; >» У. 
Теорема 1. На интегральных уравнениях Гаммерштейна вида (1) алгоритм 1 экви- 
валентен алгоритму 2. 
Доказательство. Начальные приближения у»о = 0 алгоритмов | и 2 тождествен- 
ны. Следовательно, на первой итерации (5 = 1) алгоритм 1 вычисляет: 
— функцию Д, тождественную функции / алгоритма 2, 
— многочлен А! (2), тождественный многочлену Ё\ уравнения (15), 
— многочлен 1 (3), тождественный многочлену 21 уравнения (15), 
— уравнение (4), тождественное уравнению (15), 
— многочлен уж.1 (6). 
Согласно теореме 1 [2], многочлен у„1 (6) тождественен многочлену у» (16). 
На следующих итерациях эти тождества сохраняются. 


4. Сходимость алгоритма 1 


Из теоремы 1, оптимальности операторов 5„ и (И, как аппарата аппроксимации 
в пространстве [2(а,Б, р) и Саы [3, с. 77-95] 


Пу $=[ У ] [2ьь) = Не 0,...с п | УС9 — (со +...4 си Хх”) | 12а, 
| $ зав =1, [15 [сшы = (4/ п”) ви) + О), О) <3, п>0, 
| Ч, |ааы = (2 / п) Ш(и) + О(0), О@)<1 п>0 
можно сделать следующее заключение. 
Вывод 4. Результаты исследования метода Галёркина [4, гл. 4] на достаточ- 
но широком классе уравнений Гаммерштейна (1) доказывают: 
— существование многочлена у, (6) — аппроксимации искомого решения, 
— сходимость последовательности многочленов (6) к точному решению уравнения 


у= Пти —о0о Уп › 
— Ограниченность коэффициента оптимальности алгоритма 1. 


5. Оценки погрешности алгоритма 1 


Апостериорные оценки. Теорема 3 [2] устанавливает точные и конструктив- 
ные апостериорные оценки нормы в пространстве Си» погрешности решения 
ЛИУМК (4) по алгоритму [2]. Если последовательность этих оценок для 5 = [,2,... 
сходится, то ее предел естественно принять за оценку погрешности алгоритма 1 на 
уравнении Гаммерштейна (1). 

Пример 1. Для решения уравнения (7) по алгоритму 1 мы преобразовали это 
уравнение в уравнение 


У- КГу]= КГ ехро») -у ], (17) 
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где КГ у | - линейный интегральный оператор уравнения (7). ЛИУМК (4) для уравне- 
ния (17) имеет вид 


У-Ку] = К] Ч ехр(,:-1) — Улиз-1 ] ]. (18) 
Оценка по теореме 3 [2] нормы в пространстве Суо,1] погрешности решения по алго- 
ритму [2] ЛИУМК (18) на отрезке [0,1] с параметром и имеет вид 


| И и +2(х) ||сто,и] | 7,5) |, 5 = 1,2,..., (19) 


где Ки,5) = (2 [п/2],5) — отличный от нуля коэффициент дополнительного многочле- 
на решения ЛИУМК (18) по алгоритму 1 с параметром и. Если 5 > 2, то элементы 
последовательности т, 1), (и,2),... не изменяются. пля п = 6 эта последовательность 
имеет значения 6.8 10°, 1.839 1079, 1.83857 107°,..., 1.83857 107%... Эти последо- 
вательности имеют пределы 


{ (2), (4),..., ®10) } = { 0.00028, 2.1 10°, 1.8 10%, 1.7 107°, 1.6 107? }. 


Для оператора Е У те у-К[у] ЛИУМК (18) норма функции 
И’ (х) = Г сВеБи, 2 х- 1) ] = сВеБ(я,2 х — 1) +...+ К сВеБя,2 х- 1] +... 
в пространстве Су,11 принимает следующие значения: 
| Й/> (х) ||<1о,1] — | | сНеб(2 7 2х- О] ||<1о, 1] — 1, 1= 0, в | (20) 
Следовательно, предел полученных по теореме 1.2 оценок (19) погрешности алго- 
ритма [2] на ЛИУМК (18) при 5 -> с является эффективной оценкой погрешности 
алгоритма | на уравнении Гаммерштейна (7). 

Априорные оценки. Теорема 4 [2] устанавливает точные и конструктивные 
априорные оценки нормы в пространстве Си] погрешности решения ЛИУМК (4) по 
алгоритму [2] и доказывает ограниченность коэффициента оптимальности алгоритма [2]. 
Этот коэффициент зависит только от оператора Г] у ]. Если линейная часть функции Г 
уравнения Гаммерштейна (1) мала, то оценка коэффициента оптимальности решения 
ЛИУМК (4) по алгоритму [2] достаточно точно оценивает коэффициент оптималь- 
ности решения этого уравнения по алгоритму 1. 

Пример 2. Оценка по теореме 4 [2] коэффициента оптимальности в простран- 
стве С = Сул] решения по алгоритму [2] ЛИУМК (18) на отрезке [0,1] с параметром 
п=21+1 или п=2 т имеет вид 


Си авогииа_[2], (18), С) = | № н2(% [с / а2н2( № н2 о), [0,1]. (21) 
Мы вычислили составляющие оценки (21): || И н2(х) ||сю.и=1 (20) и 
аз К И” (х), [0, |) =1+ 62» (22) 


{ дно. = { 0.06, 0.04, 0.008, 0.0036, 0.0022, 0.0013, 0.00087 1. 

Следовательно, полученные по теореме 4 [2] оценки (21) коэффициента оптималь- 
ности алгоритма [2] на ЛИУМК (18) являются эффективной оценкой коэффициента 
оптимальности алгоритма 1 на уравнении Гаммерштейна (7) и (17) в пространстве Сол]. 


6. Вычисление оценки погрешности алгоритма 1 


На основании оценок [2], мы построили алгоритмы 2 и 3 вычисления апосте- 
риорных и априорных оценок погрешности решения ЛИУМК (4) по алгоритму [2] в 
пространстве Ста]. Для ЛИУМК (18) по этим алгоритмам с параметром 4 вычисляют 
многочлен 

И’ ч(х) = ИХ, | ах) — ых) [сл < | <(и, 4) | 
и один отличный от нуля коэффициент дополнительного многочлена х(и, 4). Этот ко- 
эффициент убывает к нулю (с ростом параметра 4) 


| (и, 4) | = о(и1""), и < 0.1. (23) 
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Поведение коэффициента х(и, 4) хорошо иллюстрируют его значения 

О 101 ={210"°.3 10°. 810”. 4107, 3107”. 510. 1%. 
Следовательно, в случае ЛИУМК (18) отрезка [0,1] и параметра п: 
— вычисленная по алгоритму 2 аппроксимация 


|| 75 изза) Пет, | ®(и, $), $=1,2,... 
оценки (19) имеет погрешность 
РЕЙ изна(х) [|сто, — | 72 знач) [|столл [| <, $) [< | <(2[и/2] +2, 4) [| х(ь 5) |, 


— коэффициент Фурье — Чебышева порядка и многочлена И’,.(х) аппроксимирует ко- 
эффициент Фурье — Чебышева порядка и функции И’„(х) с погрешностью 


| ат И’, 4(х) ‚ [0,1] ) — а Их) , [0,1] ) | < | (в а) |, 
— вычисленная по алгоритму 3 аппроксимация 
| И ина.) [сон Маз ир+2( И ин, а(Х),[ 0,1] ) 
оценки (21) коэффициента оптимальности алгоритма [2] имеет погрешность 
О(| (2 [и/2 2, 4) | ), 
— коэффициент “(и, 4) удовлетворяет тождество (23). 


7. Программирование алгоритма 1 в АР5 


Структура данных на входе. 

1. Уравнение (1) определяют операторы Г, Н, функция / и многочлен . 

2. Линейные интегральные операторы Ё = Шу |, Н = ВН] у] являются суммой 
операторов вида 

Ау, [и Код О, хе [4,6] > [69,469] < [4,1 , 
где ядра, пределы интегрирования и коэффициент А являются алгебраическими 
многочленами. Эти операторы имеют вид 
А *у- шЕ ор(с_ 1, а ТК 1*у, 9+ 
+... щЕ ор(с $, 4 $, К $*у, %). 

Искомая функция у является атомом. 
Многочлены А, с 1,..., с $, Ч 1,..., 4_$ являются термами. 
Они имеют естественный для математики вид. Их аргумент х является атомом х. 
Ядра — многочлены К 1, ..., К_ $ являются термами. 
Они имеют естественный для математики вид. 
Их аргументы х, { являются атомами х, 1. 

3. Функция (7) = (ху) — нелинейность уравнения (1) — имеет естественный для 
математики вид и её аргументы х, у являются атомами х, 1. 

4. Многочлен © является термом и имеет естественный для математики вид. Его 
аргумент х является атомом х. 

5. Отрезок [а, 6] определяет список (а, Ъ) ". 

6. Параметр п алгоритма является целым числом. 

Структура данных на выходе. Многочлен у, (6) имеет числовые коэффи- 
циенты и естественный для математики вид А+... +Ё*х^п. 

АРГАМ-процедура 1. Алгебраическая спецификация п. 2 алгоритма 1. 


Е $ := $165( у=у п, Бу); ев */ 
Е $:= СВеБ_пиегро(Ё $ , ищегуа1 ‚ п); /* Ч п] $] */ 
© $:= $16 1 ц(Ну,Е $5); Ио 
ГЛОМК := (Гу+$+ (-1) * в $=0); /* ЛИУМК */ 


АРГАМ-процедура — реализация алгоритма [2]. 
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Структура выхода операторов АРГАМ-процедуры. 
На каждой итерации процедура вычисляет: 
— функцию / = Дх,ун 51) естественного для математики вида, 
— многочлен РЁ; = (| ] (2) естественного для математики вида, 
— многочлен ©; = Н] Е; | (3) естественного для математики вида, 
— ЛИУМК (4) вида, принятого в процедуре, реализующей алгоритм [2] 
А *у- шЕ ор(с1,а1,К 1*у,0-... 
+ шЕ ор(с $5,4 $, К з*у, Ю+е+ -1*е $=0, 
— решение ЛИУМК (4) процедурой, реализующей алгоритм [2], — многочлен у» с 
числовыми коэффициентами естественного для математики вида. 


8. Исследование АРГАК-процедуры 


Теорема 2. АРГАМ-процедура на каждой итерации имеет по параметру и поли- 
номиальную сложность 
(т-+1) О( сапр!№, т) + О( п’), 
т < шах { 4е2( ЦуЕР,]), 94е>( Н уЕР, |+) 1 =и- О(1), 
где О(сапрШ, т ) — сложность преобразования оператором сапрШ многочлена Р. 
Многочлен сапрЖ Р) является суммой мономов вида 


с(1) *х^]$Ъ, 1=0,...т, ]=0,...,т. 


Доказательство. АРГАМ-процедура — линейная. Следовательно, вычислительная 
сложность процедуры тождественна сумме вычислительной сложности операторов 
этой процедуры. Вычислительная сложность процедуры, реализующей алгоритм [2], 
является доминирующей. Остальные операторы процедуры имеют по параметру п 
полиномиальную сложность О(и”), $<3. 


9. Вычислительный эксперимент с процедурой 


Описание уравнения (17) на языке АРГАМ. 
ргосезз[1] := (Бу := ехр(у) + -1*у; ®:=0; 
Гу :=у+ шЕ ор(0, х, (-1 * ((х- -1) * 9) *у,® + 
ше ор(х, 1, (-1 * (х * (&+ -1))) *у,®; 
Ну := шё ор(0, х, ((х+-1) *9*у, 9 + 
ше ор(х, 1, (х * (&+ -1)) *у, ®) ; шегуа[ := (0, 1) ;;...). 
Результаты решения процедурой уравнения (17) (ип=2). 
Ё $ := 9165(у =у п, Бу) = 1; 
Е $:= Свеб пиегро!(Ёу п ‚ пщегуа|, п) = Г. 
2 $:= 3 1 (Ну, Е $) += 0.5*х^2+-0.5*х; 
ГЛОМК := (Гу+з+ -1*Е $ =0) =у+ 
ше ор(О, х, (-1 * ((х- -1) * 9) *у, О + шЕ орСх, 1, 
(-1* (х* (+ -1))) *у, © + -0.5 *х^2+0.5 *х; 
„..уУп= 0.452697 *х^2 + -0.452697 * х + -0.000294724; 
... Е_ $ := Сцеб пщегро!(Ку_п ‚ пиегуа| ‚ п) = 
-0.0248033 *х^2 + 0.0248033 *х+1; 
... уп= 0.455043 *х^2 + -0.455043 *х + -0.000280104; 
... уп= 0.455067 *х^2 + -0.455067 *х + -0.000279958. 


Вывод 5. Результаты вычислительного эксперимента -— иллюстрация эффек- 
тивности теоремы 2. 
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10. Программирование алгоритма 1 в СКА 


Построенная выше АРГАМ-процедура доказывает возможность реализации алго- 
ритма 1 в других СКА. Результаты исследования этой процедуры доказывают эффек- 
тивность такой реализации. 


11. Практическая задача 


Конечные деформации упругой струны под действием поперечной нагрузки 
хорошо моделирует краевая задача 
у"=1+а' (7), У(0)=0, У1=0. (24) 
Эта задача предложена В. Прагером [1, с. 43]. Задача (25) эквивалентна уравнению 
У=К ча? (1, Ку1= 6 а-х Ехо “+ а-дху® @& = (25) 
Задача (24) и уравнение (25) имеют единственное решение — аналитическую в 
окрестности отрезка [0,1] комплексной плоскости функцию 
у = Иа? ш( соз(а (х- 1/2) / соз(а/2) ). (26) 
Отличные от нуля коэффициенты Фурье — Чебышева функции у (26) на отрезке [0,1] 
только четные. С возрастанием параметра п = 2 Гэти коэффициенты регулярно убы- 
вают к нулю | 42 Ку, [0,1] ) | = о(4) ‚ а < 0.01 и, в случае параметра а = 0.7, при- 
нимают следующие значения: 


аку, [0,1] )}:-о ={ 0.064, -0.064 , -0.00034, -2.8 10°, -2.7 10%, 
22,1100, 52.9107”, 32107” 44610 °. 6.21077. (27) 


Следовательно, для величины наилучшего приближения функции у (26) (а = 0.7) алгеб- 
раическими многочленами порядка п в пространстве Сол] справедливо тождество (10). 

Вычислительный эксперимент с алгоритмом 1. С возрастанием параметра и 
алгоритма 1, норма погрешности решения уравнения (25) с параметром а = 0.7 по 
алгоритму 1 в пространстве Сул] убывает к нулю и удовлетворяет тождества 

{Пу- 2; |ао.и г = { 0.00036, 3.4 10°, 3.1 10%, 2.9 107%, 6.6 107?}. — (28) 

Скорость убывания тождественна скорости убывания четных коэффициентов Фурье — 
Чебышева (27) функции у (26) на отрезке [0,1] имеют место тождества 


|У— У» [“о,т / | 42 иву+2(Су › [0,1] ) |= 1+ 02 изъ { оо; = { 0.063, 0.21, 0.17, 0.08 }. 
Коэффициент оптимальности метода решения задачи С, ( тео, 145%, Свы)>1. 
Поэтому из этих тождеств и тождества (10) можно сделать следующее заключение. 


Вывод 6. Коэффициент оптимальности решения по алгоритму 1 уравнения (25) 
с параметром а = 0.7 в пространстве Су, 1] — асимптотически минимальный 


С„(аюоги т 1, 1азК (25) (а= 0.7), Со) = (1+ о(1) ) (1+ 92 21 ). 


Сравнение алгоритма 1 с методом квазилинеаризации [1]. На краевой зада- 
че (24), эквивалентной уравнению (25), согласно [1, с. 42], метод квазилинеаризации 
Из этого тождества и тождеств (27) - (29) можно сделать следующее заключение. 

Вывод 7. Коэффициент оптимальности метода (13) на краевой задаче (24) с 
параметром а = 0.7 в пространстве Сул) (эквивалентной уравнению (25)) в 
1 000 000 раз больше коэффициента оптимальности алгоритма [ на уравнении (25) 


Сто( Рио [ Ув-ло |, ак (24), Сол) > 3 10° /2.9 107? > 1 000 000. 
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Заключение 


Алгоритм 1 решает в СКА интегральные уравнения Гаммерштейна вида (1) на 
отрезке [а,6] и вычисляет алгебраический многочлен у, (6) порядка п = №. Эта ап- 
проксимация решения уравнения Гаммерштейна (1) оптимальна для символьных 
преобразований в СКА -— коэффициент оптимальности алгоритма в пространстве 
Си. ] ограничен на широком классе уравнений Гаммерштейна. 

Алгоритм 1 на уравнениях Гаммерштейна реализует идею В.К. Дзядыка. 

На основании а-метода построить эффективные численные методы решения 
функциональных уравнений. 


Дополнение 


1. Алгоритм 1' — модификация алгоритма 1. Алгоритм [2] заменен алгоритмом 
а-метода, где базис дополнительного многочлена является ортогональным базисом 
пространства Гильберта Ы\а.ь.. 

2. Для алгоритма 1' имеют место аналог теоремы 1 и аналоги применения тео- 
рем 1 —4 [2]. В них пространство [2(а,6; р) заменено пространством Ша. 

3. АРГАМ-процедура 1’ — модификация процедуры 1. Процедура, реализую- 
щая алгоритм [2], заменена реализацией алгоритма а-метода с ортогональным бази- 
сом пространства Гильберта Н/и,ь|. 
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П.М. Денисенко 


Оптимальний алгоритм для розв’язування штегральних рвнянь Гаммерштейна в системах 
комп’ютерно! алгебри 

У робот! побудовано алгебра?чний алгоритм для обчислення в системах комп’ютерно! алгебри алгебра? чного 
многочлену у, порядку п Е М№. Цей многочлен — апроксимащя розв’язку у = УС), х е [а,6в] 
1нтегрального равняння Гаммерштейна. Ця апроксимащя оптимальна в простор1 С. 


Р.М. Рет5епко 

Орта! Ао т о Зо!уше ве Наттег$ет Пиеога! ЕдиаНот$ шт Ве Сошршег Аюефга Зу$етв 
Ме сопзгасе4 Фе а|оебгатс а!2ог та Юг сотрийпе ш Фе сопариег а|оебга зузетл$ ве а|оебгалс ро]упопта! у, оЁ 
ог4ег п е №. Ты$ ро]упопта] 15$ Фе орйта! арргохипайоп Ююг фе Наттегует пиеота| едиайоп зоаЯоп 


У = У(х), хе [а,6] т Сыы: 


Статья поступила в редакцию 07.07.2008. 
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